Quantum criticality in itinerant metamagnets 
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Critical thermodynamics close to a metamagnetic quantum critical endpoint (QCEP) in a metal 
is discussed within the framework of spin-fluctuation theory. We analyze the effective potential for 
the Ising order parameter that is renormalized by spin-fluctuations and acquires a characteristic 
temperature dependence. The resulting magnetic field, H, and temperature, T, dependence of the 
magnetization, specific heat, thermal expansion, magnetostriction, susceptibility and the Griineisen 
parameter are determined, and the crossover lines in the {H, T) phase diagram are specified. We 
point out that the metamagnetic QCEP is intrinsically unstable with respect to a magnetoelastic 
coupling. 

PACS numbers: 



I. INTRODUCTION 

In metamagnets the magnetization exhibits a super- 
linear rise as function of the applied magnetic field H 
or even jumps at a first-order transition^i In the (H, T) 
plane such a line of first-order metamagnetic transi- 
tions terminates in a critical endpoint {H* , T* ) beyond 
which metamagnetism reduces to a crossover phenom- 
ena. Interestingly, an isolated quantum critical endpoint 
(QCEP) arises if the endpoint temperature T* can be 
tuned to zero by some external control parameter. 

The notion of a metamagnetic QCEP wds introduced 
in the context of the layered perovskite Sr3Ru2 07 (Ref.l2i, 
1^ to account for its peculiar metamagnetic properties. 
Here, the line of first-order metamagnetic transitions can 
be tuned to zero either using the orientation of the mag- 
netic field^ or by applying pressure In the former case, 
however, the putative QCEP is preempted by the forma- 
tion of a nematic phase,—!^ see Ref. d for a review. Sub- 
sequently, the metamagnetic crossover in the single-layer 
ruthenate Ca2-xSra;Ru04 (Ref. |^ as well as the heavy- 
fermion system CeRu2Si2 (Refs. [Tollllh was discussed 
in terms of metamagnetic quantum criticality motivated 
by the thermal expansion that is similar for the three 
compoundsji^, Note, however, that the metamagnetism 
of the latter system is in vicinity to a tricritical point 
associated with antiferromagnetic orderingj^ii^ and al- 
ternative interpretations have been also proposedJ^"— 

In fact, the presence of metamagnetic QCEPs might be 
generic to itinerant ferromagnets. It has been suggested 
that the ferromagnet transition in metals generally be- 
comes first-order upon tuning the Curie temperature to- 
wards zerojii^ The resulting tricritical point gives rise to 
surfaces of first-order metamagnetic transitions at finite 
H that terminate in a QCEP. There are indications that 
such a scenario applies to a number of ferromagnets, for 
example, ZrZn2 (Ref. IT9|) and UGe2 (Ref. 20). In a recent 
study, Aoki et al^ demonstrated that UCoAl exhibits 
a first-order ferromagnetic instability, and its metam- 
agnetic QCEP can be reached by applying a pressure 
p = l.SGPa. 

In a previous publication, we theoretically discussed 



the universal critical thermodynamic signatures expected 
close to a metamagnetic QCEP.^^ The metamagnetic in- 
stability is defined by a diverging differential susceptibil- 
ity x{T) at a finite H„i, which together with the Ising 
symmetry of the endpoint implies a number of charac- 
teristic phenomena. In the presence of a magnetoelas- 
tic coupling, one finds a diverging Griineisen parame- 
ter, a sign change of thermal expansion as a function of 
field, a{H), and a minimum in the specific- heat coeffi- 
cient ^{H). Moreover, the susceptibility, magnetostric- 
tion, and compressibility should exhibit the same meta- 
magnetic singularity and, as a result, one expect a pro- 
nounced crystal softening. 

In this work, we illustrate and compute these signa- 
tures within a specific model, i.e., the spin- fluctuation 
theory for quantum critical metamagnetism, which was 
put forward by Millis et al.—'^ extending earlier work 
by Moriyap2i and Yamadai^ We derive the detailed mag- 
netic field and temperature dependence of thermody- 
namic quantities, specify the crossover lines in the phase 
diagram, and give details for the predicted behaviour 
of thermal expansion that already appeared in a short 
publication.^^ 

The outline of the paper is as follows. In section HIl we 
introduce the model of Millis et al^^^, derive the phase 
diagram and compute the free energy density. In section 
mil we present the result for thermodynamic quantities, 
and we conclude in section ITVl with a summary and dis- 
cussion. 



II. SPIN-FLUCTUATION THEORY FOR 
ITINERANT METAMAGNETISM 

We follow Millis et al^^ and consider spin- 
fiuctuations in a metal close to a metamagnetic insta- 
bility. We concentrate on the immediate vicinity of the 
instability where the longitudinal part of the Pauli mag- 
netization M, becomes critical. The fluctuations around 
its mean value, (j) = M — (M), identifies an Ising order 
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that is governed by the potential 

(1) 



parameter 



The potential Vo('/>) has no specific symmetries and con- 
tains in principle all powers of </>. However, later we will 
use the freedom to shift the field by a constant in or- 
der to eliminate the cubic term from the renormalized 
effective potential. Close to the metamagnetic QCEP, 
the renormalized parameters in front of the linear and 
quadratic terms, i.e., h and r, are small. The linear term 
is proportional to the deviation of H from the critical 
field, h (X H — Hm- 

The potential for the Ising order parameter is pro- 
moted to a quantum critical theory by allowing for spatial 
and temporal variations. The resulting action reads22i2^ 



S 



dr / d'^x 
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with the gradient terms given by 



A(fc,w„) = k^ + 



(2) 



(3) 



where aj„ is a bosonic Matsubara frequency and the dy- 
namical exponent is z = 3. Similar to ferromagnetic 
paramagnons, the dynamics of the metamagnetic fluctu- 
ations is dominated by Landau damping due to the exci- 
tation of particle-hole pairs in the metal. In Eq.Q and 
in the following, we use dimensionless units for lengths 
and energies as well as for the Ising order parameter 0. 

The spin-fluctuation theory for metamagnetic quan- 
tum criticality rests on the assumptions that (i) the Tay- 
lor expansion of the bare potential in Eq. ([T]) exists and 
(ii) the static susceptibility is analytic in momenta such 
that indeed A(fc, 0) = fc^. This requires in particular that 
the electronic density of states is sufficiently smooth as 
the expansion coefficients of Eq. ([l} are directly related 
to the derivatives of the density of states at the Fermi 
level. For example, when the Fermi level is exactly lo- 
cated at a van-Hove singularity the effective bosonic the- 
ory ([2]) is not applicable. In addition, non-analyticities in 
a critical Fermi liquid can also been induced by interac- 
tion effects, which have been intensively investigated in 
recent years, see e.g. Refs. [26l - [28l It is established that 
for a ferromagnetic instability, i.e., a magnetic instability 
at zero magnetic field H = 0, virtual particle-hole pairs 
can mediate long-range interactions between the fiuctua- 
tions of the order parameter which in turn result in non- 
analyticities of the Landau potential that render the tran- 
sition first-orderjiS At a metamagnetic instability at fi- 
nite field Hm > 0, however, spin-flip processes are frozen 
out at low temperatures due to the finite Pauli magne- 
tization, and the Ising order parameter only couples to 
the relative density fluctuations of the spin-majority and 
-minority Fermi surfaces. It was argued in Ref. 27 that in 
this case the interaction mediated by particle-hole pairs 



is less singular, and the spin-fiuctuation theory remains 
analytic. Consequently, it is believed that the assump- 
tions (i) and (ii) hold and that the theory ^ is well 
defined. 



A. Fluctuation corrections to the effective 
potential 

The fluctuations of the field (p will dominate the low- 
temperature thermodynamics close to the metamagnetic 
QCEP. These fluctuations can be treated in a pertur- 
bative manner. It is nevertheless convenient to apply 
the method of the renormalization group (RG) to incor- 
porate these corrections in a systematic self-consistent 
fashioni^ For the problem of a QCEP, we find it conve- 
nient to apply the functional renormalization group (RG) 
and consider the RG fiow of the effective potential V it- 
self. 

The potential V at some RG scale A can be expressed 
as a line integral along an RG trajectory^Si^ 



/(Te^^',V;',(</.)e^^'), (4) 



2A^ 



where V" is the second derivative of the potential with 
respect to the field (j>. The effective dimensionality is 
D ^ d + z > 4 with the spatial dimension = 2 or rf = 3, 
and z ~ 3. The function / is given by 

(5) 

The momentum integral and the sum over bosonic Mat- 
subara frequencies is here understood to be regularized 
with some UV cutoff A. Replacing the Matsubara sum 
with an integral along the branch cut singularity and em- 
ploying a hard cutoff regularization, one obtains 



f{T,R) = --^K— de dkk'^-^ coth 
2tt dA Jo Jo 



2T 



(6) 



X arctan 



i?-|-fc2 



where = {2'^-^T:'^/'^T{d/2))~^ . 

The effective potential is obtained after incorporating 
all fluctuation corrections by integrating over the full RG 
trajectory, V = Voo- Two different types of fluctua- 
tion corrections can be distinguished corresponding to 
the separation of the / function into two parts. 



f{T,R)^foiR) + foo{T,R), 



(7) 



with 



MR) = /(0, 0) + /("^i)(0, 0)R + (0, 0)i?2, (8) 
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where /(°'") is the n*^ derivative with respect to the sec- 
ond argument. The correction arising from /q contributes 
considerably only at the initial stage of the RG flow, i.e., 
for small values of A in Eq. These corrections can be 
conveniently absorbed into a renormalization of the bare 
parameters of the potential ([T]), 



T A 



V,.en(0) = Vo(</)) 



(9) 



In the last line we disregarded an uninteresting constant. 
As announced, the freedom to shift the field allows to 
eliminate the cubic coupling term. 

On the other hand, the contribution of the fluctuations 
attributed to f^o only develops at the final stage of the 
RG trajectory, i.e., for large A and comprises the univer- 
sal fluctuation corrections. So we finally obtain for the 
effective potential 



(10) 



where we introduced the function 



J- logA/TV^ 

(11) 

The effective potential is determined by a differential 
equation ([T0| . However, it turns out that this differen- 
tial equation can be easily solved perturbatively for small 
temperatures T. The expression (ITO|) for the effective po- 
tential is the fundamental equation of this paper. 

We will need the limits of Ad for small and large ar- 
guments. For dimension d = 2, 3 the limiting behaviour 
can be summarized as follows 



Ad{x) 



-ai + a2X - 03a; 2 -I- 



as a; 



(12) 
X < 1 

a; > 1 



where the coefficients might depend logarithmically on 
X. The non-analytic term with coefficient 03 originates 
from the zero Matsubara mode. The term with coefficient 
03 will give rise to non-analytic Fermi liquid corrections. 
The explicit form of the coefficients are given in Appendix 

E 

We will limit ourselves to situations where all pa- 
rameters of the potential Vron are positive except the 
tuning parameter /i, which can have either sign. We 
will therefore not discuss first-order metamagnetic transi- 
tions, which are realized for r < 0. In the limit r = one 
obtains the QCEP. The quartic, u, quintic, U5, and higher 
order couplings are formally irrelevant in the RG sense as 
the model is above its upper critical dimension, D > 4. 
The fourth order term is however needed in order to sta- 
bilize the potential. The influence of the quintic term is 
negligible at lowest temperatures close to the QCEP and 
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/T - {h/u)M 



nonlinear 
regime 



' nonlinear 
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FIG. 1: The linear regime, uh^ , is centered around 

h — 0. The crossover lines are bent outwards from the cen- 
ter. At large temperatures the boundary follows an asymp- 
totic powerlaw in h. At a temperature T ~ (r/ii)^'''''*^^-' the 
boundaries start to repel and hit the zero temperature axis 
with an infinite slope. 



can be omitted, see Appendix lA II The most relevant 
coupling is h. It is the tuning parameter of the metam- 
agnetic transition h on H — Hm, where Hm is the critical 
magnetic field. After disregarding the quintic and higher 
order terms an Ising symmetry of the potential emerges 
with respect to the transformation h — )■ —h. This im- 
plies that the critical free energy density will be an even 
function of h, J^cr{—h) = J^cr{h), in the limit of lowest 
temperatures. 

The critical free energy density J^ci- is obtained by tak- 
ing the effective potential at its minimum value, J-a — 
V(0). In the following the solution (f> will be discussed. 



B. Classical field configuration 

The classical field configuration minimizes the effec- 
tive potential, V'{(t>) — 0. The tuning parameter h cou- 
ples linearly to the field cj) and thus acts as a force. The 
field within the potential adjusts to this force so that the 
system is in equilibrium. The restoring force which coun- 
terbalances h depends on the strength of the potential. 
In the following we distinguish a linear and a non-linear 
regime for the stiff and soft potential, respectively. 



Linear regime 2> uh^ 



In the linear regime, the effective potential is simply 
obtained by expanding (jlOp in the field (f). Retaining only 
the leading temperature correction we obtain 



R 11 



2! 



4! 



(13) 



For a short discussion of sub-leading temperature correc- 
tions see Appendix I A 11 In the linear regime, B? ^ uh"^ , 
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the quadratic dominates over the quartic term and the 
classical field configuration is given by 



(14) 



The stiffness of the potential R is defined by the equation 

R = r + uT'^°-^^/^A'a{RT-^/'') (15) 

with the first derivative of the function Ad defined in 
(|lip. Keeping only the leading temperature correction, 
R is approximately given by 



r3/2. 



quantum critical 
regime 



Fermi liquid 



R = r + u< d-5 (16) 

[ X2 T\— R > T2/3 

with positive coefficients r^, and we set z = i explicitly. 
Note that a finite temperature enhances the stiffness. For 
small arguments of A'^, the stiffness, i?, obtains a tem- 
perature dependence, T^^"^)/^. In the opposite limit, 
R r^/^, the thermal fluctuations lead to a tempera- 
ture correction of Fermi liquid type, T^. 

Taking the effective potential at the minimum, V(0), 
we obtain the free energy density. In lowest order in the 
tuning parameter h we get 

J-„-r^/M<i(i?T-2/-)-;^, for R^^uh\ (17) 

The first term is due to the Gaussian fluctuations around 
the local minimum; the second term is the energy stored 
in the potential V, and it is quadratic in h in the linear 
regime. Note that this latter contribution is temperature 
dependent due to fluctuation corrections to the poten- 
tial's stiffness, R. 



2. Non-linear regime R"^ <^ uh^ 



In the non-linear regime R^ ^ uh? the potential 
is soft and its curvature is determined by the field itself, 
V"(</)) w So we get 

V(0) « ~H + 1<P' + T^'^Aa {l<fT-^'^) . (18) 

The restoring force is now determined by the quartic 
term. In lowest order in temperature, the potential is 
minimized by 



6h 

u 



1/5 



sign(/i), R^ ^h^u 



(19) 



with mean-field exponent, 6 = 3. The resulting free en- 
ergy reads 



25/3 yl/3 



FIG. 2: The local curvature near the potential minimum, 
T- (V"(^))"/^ defines a crossover scale between a quantum 
critical and a Fermi liquid regime. The QCEP obtains in the 
limit r ^ when the quantum critical regime touches the 
zero temperature axis at /i = 0. For large \h\ 2> ^/r'-^/u the 
crossover boundary increases linearly with h. 



The term proportional to refiects the non-linear 

nature of the restoring force, and it is temperature in- 
dependent in the non-linear regime. This scaling only 
applies as long as \h\ ^ Ha where the cutoff- field hj^ is 
determined by the quintic coupling, see Appendix lA II 
Moreover, the corrections to the free energy arising from 
the temperature correction to the classical field configu- 
ration (IT9|) is sub-leading, see Appendix I A 21 

The crossover boundaries in the {h, T) plane between 
the linear and non-linear regimes are shown in Fig. [T] 
In the limit r <C uT*^-^"^'/^, when R is dominated 
by the thermal fluctuations the boundaries follow T ~ 
(/i/'u)2/('^+i). At a temperature T - {r/uf/^'^+'^'^ the 
boundaries start to repel each other and hit t he ze ro 
temperature axis with an infinite slope, h ^ ±y^r-^/u + 
0{T^). 



C. Free energy density 

The free energy is determined by the effective potential 
at the position of its minimum, Jcr = V (</>). The argu- 
ment of the function Ad in the expression for the effec- 
tive potential (ITOl) identifies a second important crossover 
scale in the phase diagram, V"(0) ^ T^/^. The local cur- 
vature is approximately given by V"(0) ~ ^(jP' with 
the stiffness R of the potential, see Eq. (IT5|) . For large 
local curvature, V"(^) ^ ^j^g g.gg energy has a 

standard Fermi liquid form, J-c-c = J^o — (T/To)^, and 
thermodynamics is conventional for a metal. In the limit 
of small curvature, V"(0) <^ T^/^, thermal fluctuations 
induce non-analytic T dependences resulting in unusual 
thermodynamic behaviour. The crossover line between 
this Fermi liquid and quantum critical regime is shown 
in Fig. H 

The intersection between the crossover lines of Fig. [1] 
and [5] define four different regimes, see Fig. [31 Note 
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— sj r' I'll \/ r'' ju h 

FIG. 3: Different regimes resulting from the two crossover 
lines of Fig.[T]and[2l In the limit r — ^ the regime IV shrinks 
to a point coinciding with the QCEP. At finite r all critical 
divergencies are cut off upon entering regime IV. 



that within regime III there is an additional crossover 
at r ~ mT^^"^^/^ depending on whether the stiffness R, 
([T5|) . is dominated by r or the thermal fluctuations. For 
dimensions d = 2, 3 the leading behaviour of the free en- 
ergy in the four regimes can be summarized as follows, 



-fi T^uh^)— -f2 



\h\ 



4/3 



,1/3 



2R 



I 

II 
III 

IV 

(21) 



where we set explicitely z = 3. The coefficients are 
all positive and depend only logarithmically on the ex- 
ternal parameters, see Appendix |5] It is important to 
realize that R is temperature dependent (jl6|) . It is in 
fact this temperature dependence which leads to some 
strong thermodynamic signatures in the vicinity of the 
metamagnetic transition. 



III. THERMODYNAMICS OF THE QUANTUM 
CRITICAL METAMAGNETIC CROSSOVER 

Thermodynamics measures the response of the free en- 
ergy upon variation of external control parameters like 
temperature, T, pressure, p, or magnetic field, H . There 
are several interesting second order derivatives of the 
free energy whose behavior is discussed in the follow- 
ing. The specific heat coefficient quantifies the change 
of the free energy upon temperature variations and is 
defined as 7 = — The thermal expansion is a 

mixed derivative with respect to temperature and pres- 
1 d T rpj^g compressibility is given 



F^l^ ~ ~F^pf ' magnetostriction is a 
mixed derivative with respect to pressure and magnetic 
field, A = -^1^ = — y QpQfi and the magnetic suscepti- 
We will also discuss the Griineisen 
ratio r = a/C and its magnetic analogue Th = ^^^/^ 



by - = -^^1 



bility is X = 



where C = 7T 
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A. Dependencies of the scaling fields 

In order to determine thermodynamics of metamag- 
netism, we have to discuss the dependence of the free 
energy (|2ip on pressure, magnetic field and temperature. 
The latter plays a special role as it appears explicitly 
in the effective theory via the dynamics of the metamag- 
netic order parameter, see Eq. ([3]), resulting in an explicit 
T-dependence of the free energy l|2ip that is responsi- 
ble for the leading low-T behavior. Besides this explicit 
T-dependence, there is also an implicit dependence via 
the scaling fields, i.e., the coefficients of the effective po- 
tential, see Eq. ([9]). They also depend on pressure and 
magnetic field. We have to consider the implicit depen- 
dence of the linear and the quadratic coefficient of the 
potential, h and r, respectively. 



h = h{H,p,T) 
r = r{H,p, T) ; 



! h{H,p) + hTT^ 
r{H,p). 



(22) 
(23) 



sure, (X y grp \p 



V OTdp ■ 



In principle, their dependence can be obtained from mi- 
croscopic calculations)^ Within our effective theory, we 
assume that the pressure and magnetic field dependence 
is sufficiently weak such that the derivatives of h and r 
with respect to p or H can be treated as a constant near 
the metamagnetic transition, e.g. dph const.. In addi- 
tion, the intrinsic temperature dependence is assumed to 
be of Fermi liquid type, T^. The temperature dependence 
of the quadratic coefficient, r, can be neglected close to 
the QCEP. However, as we will see, the dependence 
of the linear coefficient with amplitude hx leads to im- 
portant corrections although hx is formally irrelevant in 
the RG sense. 

The following discussion is divided into two parts. In 
the first part the implicit p and H dependence of only 
the most relevant parameter, h{H,p), is considered and 
/iT and the dependences of r are neglected. In this case, 
the emergent Ising symmetry of the potential jl]) will be 
directly reffected in the physical phase diagram spanned 
by temperature, T, and magnetic field, H. In a second 
step, we will consider the corrections induced by hx and 
the dependences of r{p,H). 



B. Emergent Ising symmetry in thermodynamics 

Taking into account only the p and H dependence of 
the control parameter, h{H,p) cx H — Hm{p), the critical 
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free energy is of the form, 



(24) 



Near the metamagnetic transition the control param- 
eter measures the distance to the critical field, h cx 
H — Hm{p), where the critical field is assumed to de- 
pend smoothly on pressure, dpH„i ~ const.. With these 
assumptions it is clear that the derivative of the free en- 
ergy with respect to magnetic field or pressure yield the 
same thermodynamic information. As a consequence, the 
second order derivatives of the free energy are related. 

For example, the susceptibility x is asymptotically pro- 
portional to the critical part of the compressibility k. 
In fact, their proportionality can be taken as a criterion 
for the assumption to be valid. Note that the free 
energy Tcr only accounts for the critical metamagnetic 
contributions. There might exist non-critical background 
contributions as, for example, for the compressibility. 

Another immediate consequence of the form for 
the free energy is that the emergent Ising symmetry of 
the potential ^ results in a reflection symmetry with 
respect to the metamagnetic field Hm, 

(i/ - Hra {p) , T) = J-er (i/„, (p) -H^T). (26) 

The asymptotic Ising symmetry of the problem thus be- 
comes explicit in the thermodynamic signatures. This 
symmetry implies for example that the thermal expan- 
sion as a function of H changes its sign at the critical 
field i?m, and that specific heat and susceptibility are 
even functions oi H — Hm ■ Deviations from this symme- 
try is an important indicator for additional contributions 
to the reduced form of Tcr, (EH)- 

In the following, the magnetization and the behavior 
of the second order derivatives (|25l) in the four regimes 
of Fig. [3] are discussed in detail for general dimensions 
d = 2, 3 using the expression for the free energy given in 
Section III C[ which is correct up to logarithmic correc- 
tions. Readers who are interested in these logarithmic 
corrections should consult Appendix [Bl 



1. Magnetization 

The hallmark of metamagnetism is the strong increase 
of magnetization, M, as a function of magnetic field. 
From the free energy (PTjl we obtain a mean-field be- 
haviour for the change in the magnetization, 5M = 
M — Aim, see the inset of Fig. IH^a). Close to the metam- 
agnetic field in the linear regime of Fig. [1] the magneti- 
zation increases linear as a function of H with a suscep- 
tibility determined by the stiffness R, see Eq. ([T5|). For 
larger distances, h oc H — Hm, the magnetization crosses 
over into the non- linear regime and follows SM ~ \h\^^^ 
with a mean-field exponent 6 — 3. The Ising symmetry 
is reflected in a point symmetry of the change in magne- 
tization 6M as a function of h. 
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FIG. 4: Sketch of the temperature dependence of the suscepti- 
bility —df^TcT, that is an even function of h, for (a) fields very 
close to the metamagnetic transition, \h\ < \J /u, and (b) 
larger fields such that the non-linear regime is entered for low 
T (assuming h > Q), see Fig. [S] Magnetostriction and com- 
pressibility are expected to follow the same behaviour. Inset: 
behaviour of the magnetization change 5M — M — Mm', the 
crossover lines separate the linear from the non-linear regime, 
see Fig. [T] 



2. Susceptibility, magnetostriction and compressibility 

Magnetic susceptibility, magnetostriction and com- 
pressibility asymptotically follow the same behaviour 
near the metamagnetic field and are proportional to 
—df^J^cr- They are even function of the control parameter, 
e.g. x{~h) — x{h)- In Fig.|3]a sketch of the temperature 
dependence is shown. The temperature trace in Fig.Hl^a) 



is located within the linear regime, \h\ < y^r^ /u, see 
Fig. 131 and here the susceptibility increases monotonously 
with decreasing temperature. In this regime the suscep- 
tibility is simply given by the inverse stiffness of the po- 
tential, R~^, see Eq. In Fig. HJ^b) a temperature 
trace for larger h is shown that crosses the boundary to 
the non-linear regime. At the crossover between III and 
II the behaviour changes qualitatively leading to a char- 
acteristic peak in the susceptibility as already discussed 
m The difference between the height of the peak 
and the saturation value at T = increases as the meta- 
magnetic field is approached. 

As the QCEP is approached, the derivative —d^J-a 
diverges. Note that this implies a diverging negative cor- 
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FIG. 5: Sketch of the temperature dependence of the deriva- 
tive dhdrJ'ci that is proportional to the thermal expansion a. 
It is an odd function of h, a{ — h) = —a{h), and has a charac- 
teristic maximum, (a) Behaviour very close to the metamag- 
netic transition \h\ < ^Jr^Ju and (b) for larger /i > 0. 



FIG. 6: Specific heat coefficient as a function of the tuning 
parameter /i, that is an even function of h. In panel (a) the 
behaviour is shown for temperatures within the Fermi liquid 
regime, see Fig. O Panel (b) displays the behaviour for larger 
temperatures where the quantum critical regime is entered 
with a peak at the corresponding crossover. 



rection to the compressibility! This indicates that the 
QCEP is inherently unstable in the presence of a magne- 
toelastic coupling,^!^ 



3. Thermal expansion, temperature derivative of 
magnetization 

We discuss the mixed derivative dhOx^'cr representa- 
tively in terms of the thermal expansion, a. Due to the 
Ising symmetry, the thermal expansion is an odd func- 
tion of /i, a{—h) — —a{h), and changes its sign across 
the metamagnetic transition. As detailed in Ref. [3l| a 
sign change of the thermal expansion and, as a conse- 
quence, of the Griineisen parameter, is a generic prop- 
erty of pressure-tuned quantum critical phenomena. The 
emergent Ising symmetry in the present case implies that 
the sign change is located exactly at the metamagnetic 
field H = H,n- Corrections to the positions of vanishing 
thermal expansion are discussed below in Section UlI CI 

A sketch of dhdrJ-a is shown in Fig. [5l In the upper 
panel (a), a temperature trace within the linear regime is 
shown. Here, the thermal expansion is dominated by the 
contribution of the condensate, dhdrJ-ct — —hdTR~^, 
and is thus proportional to h. First, it increases with 
lowering temperature; at the crossover T ~ (r/u)^^^''''^^^ 



within the regime III the thermal corrections cease to 
dominate the stiffness R of the effective potential re- 
sulting in a peak in a. Upon entering the Fermi liquid 
regime IV the thermal expansion vanishes linearly with 
T. In the lower panel (b) the behaviour is shown for 
\h\ > r^/u where the temperature trace crosses into 
the non-linear regime. The maximum follows here the 
boundary between III and II. Moreover, in regimes I and 
II the thermal expansion is not simply proportional to 
h anymore. The behaviour in Fig [5{b) was already dis- 
cussed in Ref. [l3 including a comparison to experiments 
on Sr3Ru2 07. 



4- Specific heat coefficient 

The specific heat is like the susceptibility an even func- 
tion of the distance to the metamagnetic transition, h. In 
Fig.[6]the specific heat coefficient, —d^J-cr, as a function 
of h is shown. In panel (a) the h trace is located within 
the Fermi liquid regime, see Fig. [31 The specific heat 
coefficient increases monotonously with lowering h and 
saturates in regime IV at a value determined by r. The 
trace in panel (b) shows the behaviour for larger temper- 
atures as one enters the quantum critical regime. Simi- 
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lar as the susceptibility and the thermal expansion, the 
specific heat coefficient as a function of h shows a char- 
acteristic peak at the boundary between the regimes II 
and III. The peak is attributed to the contribution of the 
condensate, {h? /2)d^R^^ . For small h and large temper- 
atures it yields a positive correction giving rise to a peak. 
As a function of temperature this correction changes its 
sign at T ~ (r/u)3/(''+i) and, as a consequence, for lower 
temperatures the peak is absent in the h trace. The peak 
height decreases as T = /i = is approached in contrast 
to the peak height of the susceptibility. 

As detailed in Ref [TTI. such a double peak structure 
close to the critical field is generally expected whenever 
the critical susceptibility increases with lowering temper- 
ature, dj'xiT) > 0. It follows from the Maxwell relation 
i9|f7 = d^x{T) that the curvature of the specific heat 
coefficient as a function of H coincides with that of the 
susceptibility as a function of T. Consequently, the spe- 
cific heat coefficient 7(i/) has a minimum at Hm that is 
framed by a two peaks. 



5. Gruneisen parameters 

Within the approximation (|26p , the thermal expansion 
is proportional to OtM so that the singular parts of the 
Griineisen parameter and its magnetic counterpart coin- 
cide, r (X Th cx dhdrJ-cx/iTd'^Tcr)- Due to the Ising 
symmetry, the Griineisen parameter changes sign and 
thus vanishes at the critical field Hm at any finite temper- 
ature. In regime III the behavior Th ^ h/{uT'-'^+^'^^/^) 
is obtained. In regime I the Gruneisen parameter is sim- 
ply proportional to the inverse of the tuning parameter, 
h oc H — Hm , 



drM 3-d 



1 



C 



3 H-H„ 



(27) 



as expected from general scaling arguments^^i Note, how- 
ever, that the prefactor does not obey scaling predictions 
because the QCEP for d = 2, 3 is above the upper crit- 
ical dimension. In spatial dimension d — S logarithmic 
corrections are present, 1/Th — {H — Hm) log see 
Appendix [B] 



C. Corrections to thermodynamic Ising symmetry: 
locations of entropy accumulation 

Often, the underlying Ising symmetry of a transition 
is hidden and not explicit, as for example in the usual 
liquid-gas transition, because the scaling fields are not 
simply related to the physical parameters. In the present 
case, the Ising symmetry is almost explicit close to the 
metamagnetic field near T = as the leading tempera- 
ture dependences are explicitly generated within the the- 
ory Nevertheless, there are corrections which are 
especially important whenever the critical leading con- 
tribution for some reason vanishes. This is for example 



the case for the locations of entropy accumulation and 
thus for the positions of vanishing thermal expansion as 
discussed in the following. 

As mentioned above, the reduced form of the free en- 
ergy (j26p implies that the thermal expansion changes its 
sign exactly at the critical field Hm- The locations of van- 
ishing thermal expansion or, equivalent ly, of the maxima 
of entropy'^^ in the {H, T) plane are thus a very sensitive 
probe for corrections to the reduced form In the 

following we will study the most important corrections 
which will shift the position of vanishing thermal expan- 
sion away from h — Q, i.e. H = Hm- To this end, we can 
limit ourselves to the linear regime, i.e., regime III and 
IV in the phase diagram of Fig. [31 

The most important corrections derive from the intrin- 
sic temperature dependence of h, i.e., the coefficient /it 
in the expansion ([22]) and the pressure dependence of the 
quadratic coefficient of the effective potential r = r(p). 
The latter gives rise to a finite value of the thermal ex- 
pansion a at ft, = that in regime III reads 



5a 



oc 1 



drdT 



(28) 



A temperature dependence of this form is familiar from 
second order quantum phase transitions!^ Nevertheless, 
we will not consider it further here as it is sub-leading 
compared with the one induced by /it as we explain in 
the following. 

The thermal expansion within the linear regime is 
dominated by the condensate contribution to the free en- 
ergy 



dhdT 



d h{T) 
df R 



(29) 



with h{T) ~ h + hxT^- The coefficient /it gives rise to 
an additive correction to thermal expansion. For temper- 
atures T ^ (r/M)'^/^''"'"^^ this correction simplifies to 



oa (X — 1 3 . 



(30) 



The two parameters /it and u are formally irrelevant in 
the RG sense but they combine to give a singular correc- 
tion. This contribution leads to a shift of the positions 
of vanishing thermal expansion away from H = Hm thus 
destroying the explicit, thermodynamic Ising symmetry 
(j26p . The importance of the intrinsic Fermi liquid tem- 
perature dependence for the thermal expansion near the 
critical field was noticed before in Ref. [36l 

In the presence of a finite /it, the locations of vanish- 
ing thermal expansion are shifted away from the critical 
magnetic field, h oc H — Hm = 0; these locations now 
follow 



5-d 



(31) 



for temperatures T ^ (r/u)^/(''+^). 
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IV. SUMMARY & DISCUSSION 

We discussed the critical thermodynamics of an itin- 
erant metamagnetic QCEP that is dominated by spin- 
fluctuations. These fluctuations generate a temperature 
dependence of the effective potential for the Ising or- 
der parameter, that reflects the z — 3 dynamics of the 
Landau-damped collective spin- fluctuations. The effec- 
tive potential was systematically evaluated with the help 
of the functional renormalization group, and we identi- 
fied different regimes in the phase diagram, see Fig. [S] 
We derived the H and T dependence of the critical free 
energy density analytically in the different regimes, see 
Eq. ((2T|) and Appendix [B] and discussed the resulting 
thermodynamics in section Hill 

In a previous work^ii we discussed the universal ther- 
modynamic signatures expected on general grounds close 
to a metamagnetic QCEP. The present calculations con- 
firm and illustrate these signatures within the spin- 
fluctuation model. In particular, we find a diverging 
Griineisen parameter, a sign change of the thermal ex- 
pansion a{H), and a minimum in the specific heat coef- 
ficient "y{H). Moreover, the susceptibility, magnetostric- 
tion and the compressibility share the same critical sin- 
gularities. 

This latter finding has some interesting consequences. 
The diverging critical correction to the compressibility 
implies that the metamagnetic QCEP is intrinsically un- 
stable with respect to a magnetoelastic coupling. Suffi- 
ciently close to the putative QCEP, the critical metam- 
agnetic fluctuations will strongly renormalize the elastic 
constants and destabilize the crystal. We thus expect 
that the metamagnetic QCEP will be preempted by a 
structural instability whose properties will be the sub- 
ject of a future publicationi^ 

Our results are in particular relevant for itinerant fer- 
romagnets with a tricritical point at zero fieldJ^ In such 
a case there exists a metamagnetic QCEP at some fi- 
nite Hm whose critical thermodynamics might obey the 
predictions of spin-fluctuation theory. Such ferromag- 
nets with a metamagnetic QCEP like for example UCoAl 
(Ref. 2ll) are thus promising candidates for a quantita- 
tive test of spin-fluctuation theory for critical metamag- 
netism. 
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Appendix A: Sub- leading corrections 

1. Quintic coupling 

The Ising symmetry of the potential ([SJ emerges when 
the corrections due to quintic coupling can be ne- 
glected. At non-zero temperatures, the coupling will 
generate a cubic term which in turn renormalizes the lin- 
ear coupling h. With these additional thermal corrections 
the expanded potential in the linear regime (jl3l) obtains 
the modifications, 

Vid^) = V(0) - + + ^ ^ _ ^ (Al) 

where 

-H = -h + C/3T(-°-2)/^^^(^y-2A-) 

c/3 = M5r(^-2)/^^;^(i?r-2/-) (A3) 

with the same function A'^ as in the expression for the 
stiffness R, (fTS]) . 

In regime III, see Fig. [31 the resulting temperature cor- 
rection to the tuning parameter, h, is of order 0{u5T^^^) 
in d 3 and ©(usT^ log^ T) in d ^ 2. The correction is 
thus smaller {d = 3) or of the same order (disregarding 
the logarithmic enhancement in c? = 2) as the intrinsic 
Fermi liquid temperature dependence that the parameter 
h possesses anyway. These corrections to h are thus not 
considered further. (The correction to thermodynamics 
due to an intrinsic depedence of h are discussed in 
Section IIII CI ) Moreover, it can be easily checked that 
the presence of the cubic term leads only to sub- 
leading corrections to the classical field configuration 0, 
see Eq. p^ . 

In the non-linear regime, on the other hand, the quin- 
tic coupling can only be neglected at sufficiently small 
h. If the tuning h reaches a threshold \h\ > h\ with 
/lA ~ u*/ul, corrections to the scaling of Eq. (fTOl) be- 
come important and Ising symmetry is lost. 

2. Corrections to (p in the non-linear regime 

In the derivation of the free energy in the non-linear 
regime ([20]) the zero temperature expression of the clas- 
sical field configuration 0, Eq. HH), was used. Here we 
compute the corrections arising from the induced tem- 
perature correction to the classical field. Introducing the 
deviation 5<j) ^ 0—0, the effective potential separates 
into two parts V(0) = V(0) + 5V{5(f)) with the correction 

5V{5(t))^^5^^ - f5,p. (A4) 
The temperature force is given by / = F{T) — F{{)) where 
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The temperature dependence of F possesses similar lim- 
its as the one of R in Eq. ([T5)) . By increasing the temper- 
ature one applies an effective force / onto the field that 
modifies the position of the minimum. The resulting cor- 
rection to the free energy is given by STcr = —P I 
This contribution is always sub-leading on the presented 
level of accuracy. 



Appendix B: Asymptotic expansions 

The asymptotic expansion of the function Ad is pre- 
sented, see Eq. ([TT|) . Moreover, the numerical coefficients 
for the expansion of the potential's stiffness, i?, Eq. (|15l) . 
and the free energy (PT|) are given. 



1. Dimension d = 3 



In spatial dimension d — the function Ad has the limiting behaviour 



^3(2;) 



A3 

-bi log — - bo + b2a; - ba.TJ 



for x ^ 1 



bix^ log 



A3 



boa;3 - b2 log 



A3 



bs a; log x 



for 



T2/3 



(Bl) 



> a; > 1. 



There is a reminiscent dependence on the cutoff A. The coefficients bo and bo can be absorbed in the non-universal 
logarithmic dependence of the leading behavior, and both are omitted in the following. The other coefficients read 



bl = 



1 

36^' 



b2 = 



1 

36^' 



bs = — . 

^ 60 



The leading temperature correction to the stiffness is given by 

b2 u T^/^ for R < T2/3 
for R > T2/3. 



R = 



The limiting form of the free energy in the four different regimes of Fig. |3] reads 

V2A3 . 2 |/i|2/3^1/3 (3|/i|)4/3 

"37:^^ log 



-b2r^ log 



T2/3 25/3^1/3 



-biT^ log 
A3 



A3 _ _^ 
Y ~ 2R 



rp4 



&27^'l°gl37J + ^3-jlO; 



T2/3 2i? 



I 

II 
III 

IV 



(B2) 



(B3) 



(B4) 



2. Dimension d = 2 



In spatial dimension d — 2 the asymptotic form of the function Ad is 

Af\ ) -ci + C3a;logi -f ... for a; < 1 

^2(2;) = S j_ (B5) 
Ci — C2 a:~2 -I- csx"'^ + . . . for a: 3> 1 

The function is universal in the limit A/r^/3 ^ qq. The coefficients are 

ci = -^r(5/3)C(5/3), C3 = ^, li^^, C2^^, C3 = ^C(3). (B6) 
An 8n dOn 24 An 

The term with coefficient C3 derives from the zero Matsubara mode whose contribution is logarithmically enhanced 
in d = 2. This logarithmic dependence is reflected in the temperature dependence of the stiffness 

, , .C3.Tlog2^^ for i?«TV3 
^ - • h~^2U^ for R:^T^/\ 
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For T r/u m the quantum critical regime, the temperature determines the stiffness so that asymptoticahy R 
z^uT\og{l/{uT'^/^)). The free energy simphfies in the four regimes of Fig. [3]to 



-C2 



21/2 J.2 



C3: 



T3 



(31/^1)4/3 



61/3 (if/i2) 1/6 ' ^64/3 (w/i2)2/3 25/3^1/3 
, y5/3 , 6^/3 /3 (31/^1)4/3 

_CiT5/3 _ ZL 

2i? 



-C2 



A/2 



2R 



I 

II 
III 
IV 



r 



(B8) 



The terms of the free energy with coefficients b3 and omitted in the body of the paper. 
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